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The energies of different angular momentum states of a heavy- light meson were measured on a lattice in [Q. 
We have now repeated this study using several different lattices, quenched and unquenched, that have different 
physical lattice sizes, clover coefficients and quark-gluon couplings. The heavy quark is taken to be infinitely heavy, 
whereas the light quark mass is approximately that of the strange quark. By interpolating and extrapolating in 
the light quark mass we can thus compare the lattice results with B and B„ meson experiments. Most interesting 
is the lowest P-wave B 3 state, since it is possible that it lies below the BK threshold and hence is very narrow. 
Unfortunately, there are no experimental results on P-wave B or B s mesons available at present. 

In addition to the energy spectrum, we measured earlier also vector (charge) and scalar (matter) radial distri- 
butions of the light quark in the S-wave states of a heavy- light meson on a lattice Now we are extending the 
study of radial distributions to P-wave states. 



1. Energies 

The basic quantity for evaluating the energies 
of heavy-light mesons is the 2-point correlation 
function — see Fig. ^ a ) ■ It is defined as 

C 2 (T) = {P t TG q (x,t + T,t) 

■P t+T ^U Q (yL,t,t + T)), (1) 

where (x,t,T) is the heavy (infinite mass)- 
quark propagator and G g (x, t+T, t) the light anti- 
quark propagator. P t is a linear combination of 
products of gauge links at time t along paths P 
and r defines the spin structure of the operator. 
The (...) means the average over the whole lat- 
tice. The energies are then extracted by fitting 
the Ci with a sum of exponentials, 

»„,„ 

C 2 (T) « £ c t e~^ T , (2) 

i=l 

where 7V max = 2 - 4, T < 10. 
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Figure 1. a) Two-point correlation function Ci\ 
b) Three-point correlation function C3. 
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Table 1 

Lattice parameters. 





a 


Csw 


K 


Q3 


5.7 


1.57 


0.14077 


DF2 


5.2 


1.76 


0.1395 


DF3 


5.2 


2.0171 


0.1350 




a [fm] 


m q /m s 


r m v 


Q3 


0.179(9) 


0.83 


1.555(6) 


DF2 


0.152(8) 


1.28 


1.94(3) 


DF3 


0.110(6) 


1.12 


1.93(3) 
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Figure 2. The energy spectrum for different an- 
gular momentum states. Here L+(— ) means that 
the light quark spin couples to angular momen- 
tum L giving the total j = L ± 1/2. 2S is the first 
radially excited L = state. The spectrum was 
measured using different lattices: Q3 (a 16 3 x 24 
quenched lattice, [T]) and DF2, DF3 (16 3 x 24 and 
16 3 x 32 lattices respectively, dynamical fermions, 
PJ). Here r$ = 0.525(25) fm. Lattice parameters 
are given in Table ^ Note that the spectrum 
shows an approximately linear rise in excitation 
energy with L. 
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Figure 3. Interpolation from infinite mass to 
mQ = nib. The results at vtlq — m c are from ex- 
periments. Note that the lowest B s P-wave state 
lies below the BK threshold and should be very 



narrow. 




From some theoretical considerations 0] it is 
expected that, for higher angular momentum 
states, the multiplets should be inverted com- 
pared with the Coulomb spectrum, i.e. L— 
should lie higher than L+ (see caption of Fig. [5] 
for notation) . Experimentally this inversion is not 
seen for P-waves, and now the lattice measure- 
ments show that there is no inversion in the D- 
wave states either. In fact, the D+ and D— states 
seem to be nearly degenerate, i.e. the spin-orbit 
splitting is very small (see Fig. |5J). 

2. Radial distributions 

For evaluating the radial distributions of the 
light quark a three-point correlation function is 
needed — see Fig.^b). It is defined as 



C 3 (R,T) = (r^U Q TG ql eG q2 ). 



(3) 



We have now two light quark propagators, G q \ 
and G q 2, and a probe 0(i?) at distance R from 
the static quark. We have used two probes: 74 for 
the vector (charge) and 1 for the scalar (matter) 



distribution. The radial distributions, x' 
are then extracted by fitting the C3 with 



C 3 (R,T) 



L 

i,j=l 



(R)e 



(4) 



— see Figs.0J|51 The m^s and c^'s are from Eq.[21 
The lattice measurements, i.e. the energy spec- 
trum and the radial distributions, can be used 
to test potential models. There are several ad- 
vantages for model makers: First of all, since 
the heavy quark is infinitely heavy, we have es- 
sentially a one-body problem. Secondly, on the 
lattice we know what we put in — one heavy 
quark and one light anti-quark — which makes 
the lattice system much more simple than the real 
world. We can, for example, try to use the Dirac 
equation to interpret the S- and P-wave distri- 
butions. Reasonable parameters can give a sur- 
prisingly good qualitative fit to the distributions 
measured on the lattice — see Fig. El 

3. Main conclusions 

• There should be several narrow B s states, 
e.g. + that lies below the BK threshold. 
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Figure 4. Ground state vector and scalar radial 
distributions and discretized exponential fits. 




Figure 5. Radial distributions containing the first 
radially excited state (Ref. for the S-wave). 




Figure 6. Dirac equation numerical solutions for a 
potential V — ~e/R + b- R compared with lattice 
results. Here e = 0.6 • He, b = (500 MeV) 2 and 
m — fOO MeV. R is in lattice units (a ps O.f 5 fm). 
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• The spin-orbit splitting is very small. 

• The radial distributions of S and P— states 
can be qualitatively understood by using a 
Dirac equation model. 
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